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Abstract 

Leptonic decays of 5-mesons are theoretically very clean probes for testing the 
Standard Model (SM) and possible physics beyond it. Amongst the various 
leptonic decays of the 5-meson, the pure dileptonic decay B — ► £ + £~ is very 
important as this mode is helicity suppressed in the SM but can be substantially 
enhanced in some of the models beyond the SM, such as supersymmetric (SUSY) 
theories and the two Higgs doublet model (2HDM). Although the purely dilep- 
tonic decay mode is helicity suppressed in the SM its associated mode B — > £ + £~7 
does not have the same suppression, due to the presence of the 7 in the final state. 
In this paper we will also analyse the effects of enhanced Z-penguins on these 
two decay modes. 



*src@physics.du.ac.in 
talanc@yukawa.kyoto-u.ac.jp 
t naveen@physics . du . ac . in 
§ j oshi@physics . unimelb .edu.au 



1 Introduction 



Experimental data on the decay of the 5-meson into two pseudo-scalar mesons seems to 
indicate a "puzzle" , in the sense that it points to a sizeable deviation from standard theories, 
as pointed out in Buras et aJ.[l]. To understand the origins of this "puzzle" consider the decay 
B — > 7T7T, where this decay is reasonably well described within the theoretical framework of 
an effective Hamiltonian [2]. However, extensions of these results using SU(3) symmetries 
for decays B — > irK show considerable disagreement with data [1-3]. A possible resolution 
of the "B — > irK" puzzle by Buras et ai.[l] has aroused considerable interest [1-3]. The 
explanation is to attribute this deviation to an enhanced Z-penguin diagram in addition to 
this diagram having an additional large phase. This explanation has aroused much interest 
as it could be the first indication of physics beyond the standard model (SM), where it has 
been usually accepted that the values of the CKM matrix cannot reproduce such a large 
phase in the Z-penguin diagram. 

The idea of a strongly enhanced Z-penguin contribution is not a new one. It was first 
carried out by Colangelo and Isidori [4] in relation to the K — > -nvv and K — > tt££ decays. 
Note that Buras and Silvestrini [5] placed constraints on these contributions in a general 
class of super symmetric (SUSY) models. The possibility of nonstandard Z-couplings in the 
context of b — > sit transitions was studied in [6]. The possible implications of this were 
studied in some of our earlier works [7]. However, [1] were the first to relate the possible 
enhancement of the Z-penguins to the nonleptonic decay modes involved in the ll B — > nK" 
puzzle and were able to obtain definitive phenomenological values of the magnitude and 
phase of the Z-penguin that is consistent with the B — > tttt, irK data. The estimates of the 
magnitude and phase of the Z-penguin required to fit the tttt and tcK data, which have been 
made in reference [2], are purely on a phenomenological basis. On the theoretical side, such 
an enhancement of the penguin diagrams can be accommodated within SUSY extensions of 
the SM, the minimal supersymmetric standard model (MSSM) in particular. The flavour 
rotation of the squarks is different in such theories from the corresponding flavour rotation 
of the quarks, and this mismatch becomes the source of an additional phase in flavour 
changing amplitudes. Any attempts to fit the evaluated value of the Z-penguin diagrams 
with theory, however, is hopeless, since the parameters involved in estimating the resultant 
phase, which are essentially the off-diagonal elements of the squark mass matrix, are not 
known [8]. Irrespective of this, the occurrence of a phase (other than the CKM phase) in 
the Z-penguin is a signal for a new source of CP violation, and has wider implications. 

The basic vertex involved in the analysis of Buras et aJ.[l] is the bsZ vertex. Processes 
which feature this vertex have been widely studied in the context of many other semi-leptonic 
and hadronic processes [1-3,7]. In this work we will address the implications of the large 
bsZ phenomenological phase to the processes B — > £~£ + and B — > £~£ +/ -f. 

Dileptonic decays of B s are very special rare leptonic decays; the non-radiative mode 
B s — > £ + £~is helicity suppressed in the SM whereas the radiative mode B s — > £ + £~^&oes 
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not suffer from such a suppression in the SM. Thus, the later can be comparable to the 
former despite an extra factor of a (the electromagnetic coupling constant). Recently the 
CDF and DO collaborations [9] have reported a very useful bound on the pure dileptonic 
mode and in the near future we may expect data on both the radiative and non-radiative 
dileptonic modes. This motivates us to look afresh for the signatures of new physics in these 
two decays. 

Within the SM, the pure dileptonic decay is dominated by the Z-penguin and box dia- 
grams, both of which are helicity suppressed. It was first recognized by two of us [10] that 
the situation changes in SUSY theories if the parameter tanf3 is large, and this has been the 
subject of investigations in a large number of later works [11, 12]. Similar sorts of enhance- 
ments are also possible in two Higgs doublet models, as has been emphasized in many earlier 
works [13, 14]. It was also pointed out by Handoko et aJ.[15] that possible measurements of 
the longitudinal polarizations in the pure dileptonic case will be a useful signature of new 
physics when SUSY effects with large tanf3 are involved. This polarization is sensitive to 
the phase of the effective bsZ vertex and hence its measurement can also provide us with 
insights into new physics. 

When a photon is emitted in addition to the lepton pair in a purely dileptonic decay, 
no helicity suppression exists. This process, B — > £ + £~'~f, is thus of special interest among 
the rare -B-decays. There have been several investigations of the sensitivity of possible 
measurables, like forward-backward (FB) asymmetries and polarization asymmetries on new 
physics [16-21]. Although this process is, theoretically, somewhat less clean as compared to 
the pure dileptonic mode, in so far as it requires a knowledge of the B s — > 7 form factors, this 
process has a greater number of observables than the pure dileptonic decay. An enhanced 
bsZ coupling with a phase will affect all these parameters, and motivates us in to calculating 
the various observables in the radiative dileptonic mode. There is also a second reason for 
reinvestigating this decay. In our earlier work we used form factors which did not meet 
certain formal constraints, as established by Kriiger and Melikhov [16], who suggested in 
their paper a recalculation of our results with a new set of form factors. We have also done 
this here with the inclusion of a non-standard bsZ coupling, as stated above. 

In the present work we will examine the pure dileptonic decay mode (B s — > £ + £~) and 
the radiative mode (B s — > £ + £~^y) in two different models, namely the two Higgs doublet 
model (2HDM) and the minimal Supergravity (mSUGRA) model. As such, this paper shall 
be organized along the following lines. In section 2 we shall present the effective Hamiltonian 
for the quark level transition b — > s£ + £~, before specialising to the purely dileptonic decay 
(B — > £~£ + ) and deriving its matrix element and other observables associated with it. The 
final subsection of this section of the paper shall analyze the radiative process (B — > £~£ +/ -f), 
deriving expressions for its matrix element, branching ratio and FB asymmetry. The paper 
shall be concluded in section 3 with our numerical results, discussions and conclusions. 
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2 The Effective Hamiltonian 



The processes B s — > i + i and B s — > £ + £ 7are based on the quark process b — > s£ + £ which 
can be described by the effective Hamiltonian; 



H 



eff 



V27T 



v tb v* 



2C 7 ^ (sia^P R b) i-fl + C 9 e// (s^P L b) t-fi 



+C 10 (s^P L b) Ind + C Ql (sP R b) It + C Q2 (sP R b) £ l5 i 



(1) 



where g M = (p + and Pl/r = |(1 =F 7s)- The values of the various Wilson coefficients 

are given in [13, 22] for the 2HDM model and [12] in the SUSY extensions of the SM. Cf f 
has the value; 

Cf = C 9 + Y(q 2 ), 

where Y(q 2 ) has contributions from one loop matrix elements of four quark operators. This, 
along with short distance contributions, also has long-distance contributions due the cc 
resonant states. These contributions we have incorporated using the prescription given in 
Kriiger and Sehgal [23] . Using the above definition of the effective Hamiltonian we will derive 
expressions for the observables we want to study. 



2.1 The Dileptonic Decay (B s -> £ + £~) 

To evaluate B s — > £ + £~we need hadronic matrix elements of various currents between B s and 
the vacuum. These hadronic matrix elements, using the PCAC ansatz, can be formulated 

as; 

{®\qi^b\B(p B )) = -if B p BfX , 
ifB~m R 

(0\q l5 b\B(p B )) =r 



m b 



(0\qa^P Rl5 b\B(p B )) = 0. 



(2) 



From the effective Hamiltonian given in eqn.(2), and using the above definition for the 
form factors, we can write the matrix element as; 



M = if Bq ^V tb VZ 



m, 



2m e C 10 + -?-C Q2 )£ l5 £ + 
m b J 



m B n 

m b 



(3) 



From the above matrix element we can evaluate the branching ratio as; 
B(B S -> £+£-] 



GW 



iy" r i ,2 ,j /-, 4m^ 

V tbV ts \ JB m B\l 1 - "^2~ 



64vr 3 



ml 



m 



2m e C 10 + ^C Q 



In the pure dileptonic (B s — > £ + £~) decay we have only one momenta available in the 
final state. Hence, we can have only one polarization asymmetry namely the longitudinal 
one. This polarization asymmetry provides a direct means for measuring the scalar and 
pseudoscalar type interactions, which are induced in almost all variants of the 2HDM and 
SUSY models. 

Defining the longitudinal polarization asymmetry of the final state leptons as [15]; 



A ± 



[T(sj, sj) + T( T sJ, ±4)] - [F(±sj, Tsf) + r(-^ 



LP 



[T(±si,Tst)+n-si,-sl)} 



(5) 



Using the matrix elment given in eqn.(3) we can derive the longitudinal polarization asym- 
metry of t~ as; 



•Alp — 



2J1 



Ami 



Re 



m B (J* 
r«5 Ql 



2m e C 10 - 



2m e C 10 + ^C Q . 



+ 



\ ~™B ) 



2 ' 



(6) 



where A£p = A[ P . Note that within the SM Alp(B s — > is zero. Therefore this 

asymmetry will be a very useful probe of the new operators in the effective Hamiltonian, in 
particular to the non CKM-phases in the various Wilson coefficients. 



2.2 The Radiative Decay (B s -> £+£->y) 

The radiative decay B s — > £ + £~7can be obtained from the quark level transition b — > s£ + £~ , 
as given by the effective Hamiltonian in eqn.(2). To obtain B s — > £ + £~7we have to attach 
a photon line to any of the charged internal and external lines of b — > s£ + £~ . As pointed 
out in [17,24] contributions coming from the photon attached to charged internal lines can 
be neglected. There are two ways to attach the photon line to the external lines; firstly to 
the external hadronic lines, and secondly to the lepton lines. These shall now be analyzed 
separately. 

The task of attaching the photon to the hadronic line can be accomplished by using the 
following form factors, as defined in Kriiger and Melikov [16]; 

{l{k)\si^b\B{p B )) = ie [e* (p B • k) - (p B • e*) fcj — , (7) 

m B 

( 7 (A;)|^6|S(p B )) = ee* a e^ apa p B k^, (8) 

m B 

(l(k)\sa^ l5 b\B(p B )) = e [el (p B ■ k) - (p B ■ e*) k,} F TA , (9) 

{^sa^Bips)) = iee^e^pPgk'FTv, (10) 

where e M and k^ are the polarization and four momenta of the photon. The definition of the 
form factors is given in Appendix A. We can express the matrix element for this first decay 
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process as; 



Mt = t= V tb V* 



2tt 



hn) {At* a e, ap(T p B k° + iB { e ; (p B ■ k) - (p B ■ e*)}) 
+ {Ce* a e^ pa p B k° + iD {el (p B ■ k) - (p B ■ e*)}) 



, (11) 



in writing the above we have used e = 2\pi\ct. And where we have used the following 
appropriately defined constants; 

A = 2Ci—F T y + C 9 — — , 
q 2 m B 

B = —2C 7 —^-F TA — C 9 —^-, 
q z m B 



C = C 



10" 



D = -C 



'v_ 

m B 

F A 



10" 



m B 



(12) 



The second process we need to consider is where the photon is attached to the lepton 
lines. This will mean making the following substitution to the lepton operators in our initial 
Hamiltonian; 



1 



_+ }t — VTii 



0\e + ie£[0 



1 



and using the form factors given in eqn.(2). 

We can therefore express the matrix element for this second process as; 



M 2 = ^V tb V* 



2tt 



[2m e C 10 + ^Cc 



-C Ql < 2m E 



m b 



+ 



1 



2p + ■ k 2p_ ■ k 



tdt + t 



m b " i \2p+ • k 2p_ ■ k 
The total matrix element would be; 

M = M! + M 2 , 
hence squared matrix element can be written as; 

\M\ 2 = \M 1 \ 2 + \M 2 \ 2 + 2Re(M 1 M* 2 ), 

where; 

a 3 ' 2 G F 



2p + ■ k 2p_ ■ k 



(13) 



(14) 



(15) 



Mi 



{\A\ 2 + \B\ 2 ) {2m 2 ((k ■ p + f + (k ■ p-f + (k ■ P+ )(k ■ p_)) 
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+(k ■ P-) 2 P- -p+ + {k- P+) 2 P- ■ P+} 

+ (|C| 2 + \D\ 2 ) {-2m 2 (k ■ p_)(k ■ p + ) +P-- P+ {(k ■ p + f + (k • p_) 2 ) } 



+Re(A*D + B*C)q 2 { (k ■ p + ) 2 + (k ■ p_) 2 } 



(16) 



\M 2 \ 2 



a 



2tt 



2 f2 
B 



Cl0 + _™5_C<; 

2mtm b 



2 + 



/c • p + + 2p_ ■ p + — m 



k ■ p- 



k ■ p- + 2p_ • p + - m| <f{p~-p+) m}(m} + k ■ p + + p + ■ p_) 



k-p+ 

m 2 e (mj + k ■ p_ + p + ■ p_) 



(k-p_)(k-p + ) 



(k-p+y 



m 



B 



2m e m b 



(k-p-) 2 

k ■ p + + 2p_ ■ p + — m 
k • p_ 



k ■ p- + 2p_ ■ p + - m| (p_ -p+ - m|)(p_ -p + ) m|(m | - k ■ p + - p + ■ p_) 



k-p + (k-p-)(k-p+) (k-p-) 2 

m 2 {m 2 — k • p_ — p + ■ p_ 



(17) 



and 

2 Re(MiM* 2 ) = 



a 



3 / 2 G, 



2tt 
(q ■ kf 



-v tb v; s 



16m 2 / f 



Re C 10 + 



m 



2mgm b 



C Q2 U* !> x 



(P+ ■ k)(p- ■ k) 



+ Re I Cio + 



mi 



+Re 



m. 



-C 0l £?* 



2m e m b 
1 



(g • k) 2 (p^ • k — p + • k) 
(P+ ' • fc ) 

- (p- " A:) 3 



+Re 



2mgm b ^ 1 J J (p+ ■ • fc) 

+ • k) 2 (2m 2 — 3p + ■ k) — (p + ■ k)(3p- ■ k + 4p„ ■ p + )(p_ • k) 

+ (p+ ■ k) 2 (2m 2 -p + ■ k) 

*\ (q-k) 2 (p_-k-p + -k) 



m, 



CqAC 



2mtm b 

The differential decay rate is then; 



dT m\ 



(p+ • fc)(p_ • fc) 



(18) 



ds 2 9 7T 3 



3/2 G; 



A, 



(19) 



with; 



A = im|(l-^1-^ 



{s + 2m 2 ) {|A| 2 + \B\ 2 } + (s- Am 2 ) {\C\ 2 + |L>| 2 } 



6 



-64 



Pb 



m. 



wo + 



2mtrrib 



2 „ o 

mi 



'1 - s 



(4m| - s - 1) log(Z) + 2s\ 1 



4m/ 



+16 



Qi 



m 6 



* Wm^-Wl ""' 



(I6m^ - 4mj(l + 2s) + s 2 + l) 



MZ)} - 32m 2 (l - s) 2 f B Re {a* (c 10 + ||^) } 



J6 ™« R e {iTC Ql } \ 4(m 2 + I - 1) log(Z) - 2 Jl - }> (20) 

m;, I V s 1 



where A (a, 6, c) = a 2 + 6 2 + c 2 — 2a6 — 2ac — 26c and Z 



1 + 



i- wi-lM 



Next we compute the FB asymmetry associated with the final state lepton. The definition 
of the FB asymmetry is; 



A FB (s) 



f 1 rf C0 ^/^ - f° dcos9-^„ 

JU ascosO Jl ascosa 



Jo acuhV dscose ^ Ji acuhV dscose 



(21) 



where 6* is the angle between £ and 7 in the dileptonic cm frame. Using this definition of 
FB asymmetry we arrive at the expression for B — > £ + £~^ as; 



A FB = - 



2m|s(l - s) 3 ( 1 



4m i 



+32 f B mj(l - syLog 



Ami 



Re (A*D + B*C) 



Re 



m 



D* ( Cio + - — — Cq 2 



+ 16_^m/(l - s) 2 Lo# 

m 6 V s 



' jA7r ' U fsRe[C*C Ql ] 



(22) 



3 Numerical results and discussion 



In this section we will discuss the results of the numerical analysis we have carried out in 
both the pure dileptonic and radiative decays. For this purpose we have divided this section 
into two subsections devoted to the numerical analysis of the two different decay modes. 
The input parameters used for our numerical estimates are given in Appendix B. 

As we have earlier stated, we are interested in estimating the effects of a large bsZ 
coupling in dileptonic decays. This sort of enhancement having been proposed by Buras et 
aJ.[2] in order to solve the B — > rnr and B — > ixK puzzle, where it was proposed that the 
effective bsZ vertex not only be enhanced in magnitude (by more than twice as much as 
compared to the SM value) but also receive a large phase, making its effective coupling to be 
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predominantly imaginary. Their fitting effectively makes Cio complex, and can be defined 

as; 

Cio = -{2.2/ sin 2 9 w )e i6Y , 6 Y = -(100 ± 12)°. (23) 

We will be using this value for Cio in our analysis with the central value of Qy. 

We have carried out our analysis with the SM using the enhanced value of Cio above. We 
have also used this enhanced value in our numerical work for the 2HDM and SUSY models. 
Note that both the 2HDM and SUSY models have the universal feature, in connection to 
the rare decays, that if we are in a large tanfi 1 region of their parameter space then we have 
to introduce " 'nonstandard 1 operators, namely the scalar and pseudo-scalar operators. In 
our case this means the introduction of operators corresponding to the Wilsons operators 
C Ql and Cq 2 in eqn.(2). 

For the case of the 2HDM we have worked in the type-II of this model. The Higgs sector 
of the type-II 2HDM has almost similar structure and Yukawa couplings to the MSSM. The 
details of this model have been given in many earlier works [13,22]. In the 2HDM type-II 
we have used the values of the Wilsons as given in [13,22]. The experimental observation 
of b — > 57 provides us with a constraint on the charged Higgs mass of this model, that is 
m H ± > 350 GeV [25]. 

Supersymmetry (SUSY) is one of the most favoured extensions of the SM. But as is 
well known, even the minimal SUSY extension of the SM, known as MSSM, has a very 
large number of parameters, making it difficult to do any phenomenology. However, the 
vast parameter space of the MSSM can be reduced to a manageable level by assuming some 
unification of the parameters of the MSSM. There are many such models generically named as 
constrained MSSM (CMSSM) available. Among those one of the very widely used models is 
the minimal Supergravity (mSUGRA) model. This model has only five parameters, namely; 
m (the unified mass of all the scalars), M (the unified mass of all gauginos), A (the unified 
trilinear couplings), tanf3 (the ratio of vev's of the Higgs doublets) and the sign of //. We 
have performed our analysis of the dileptonic decays with the mSUGRA model also. 

3.1 Pure dileptonic decay (B s — > £ + £ ) 

The expressions for the branching ratio and the longitudinal polarization asymmetry for the 
final state lepton are given in eqns.(4) and (6) respectively. Note that Buras et ai.[2] have 
already noted the enhancement in the branching ratio of B s — > fi + fi~by using the new bsZ 
vertex. 

The enhancement of the branching ratio of the pure dileptonic decay mode (B s — > /U + /z~) 
in the presence of the new set of scalar and pseudo-scalar operators has been emphasized 
earlier [10-14]. Handoko et aJ.[15] argued that these operators in the dileptonic decay (B s — > 
also provides a non-zero value of the longitudinal polarization asymmetry for the final 

1 tan[3 is the ratio of the vev's of the two Higgs doublets which are present 
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state lepton. This can also be used as a direct measurement of the physics of scalar and 
pseudo-scalar interactions. They also presented the correlation between the polarization 
asymmetry {Alp) and the branching ratio of this mode. At this point we would like to note 
that this correlation will be modified substantially if we have a large complex phase in the 
electroweak penguins as suggested by Buras et ai.[2]. Therefore the polarization asymmetry 
in the dileptonic mode can not only serve as a test of the new set of operators, but it can also 
provide us with a useful insight in to the nature of the electroweak penguins. To substantiate 
this point we have presented two plots of Alp in figure 1. The first plot being Alp as a 
function of the charged Higgs boson mass (m^) for various values of tan/3. In the second 
plot we have given the same plot in the mSUGRA model for various values of the unified 
scalar mass (m). For both the plots we have taken the value of C\q as given in eqn.(23). As 
we can see from these plots there is a marked change in the behaviour of Alp as compared 
to the case where the bsZ vertex is of the SM type. For the set of input parameters we have 
chosen we calculate the SM branching ratio as Br(B s —>■ = 3.41 x 1CT 9 . Using the 

value of Cio as given by eqn.(23) we get Br(B s — > fi + u~) = 17.21 x 1CT 9 . Such that the 
ratio of the branching ratios is; 

Br(B 8 ^n+fj,-) ^ 5 
Br(B s -> sm 




Figure 1: B s -> fi + fi~: 2HDM (left) and SUSY (right). The shaded region in the 2HDM 
figure (left) indicates the region ruled out by the b — *> 57 observations. 



3.2 Radiative decay (B s £ + £~j) 

The radiative decay in the SM does not have the helicity suppression of the pure dileptonic 
decay. This feature makes this decay very useful. Aside from this, in the radiative mode 
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one can measure a large number of observables, such as the FB asymmetry [16-19], and the 
single and double polarization asymmetries [21] associated with the final state lepton pair. 
For calculations of the observables in the radiative mode we require the definition of the form 
factors. In the literature there have been many definitions of the form factors required for 
the B — > 7 transition. The FB asymmetry for B s — > £ + £~7shows a very strong dependence 
on the specific form of the form factors being used. Where the literature provides us with 
form factors based on QCD sum rules [17], the quark model [18] etc. However, as argued 
by Kriiger and Melikov these form factors contradict some of the basic properties. They 
further gave another parameterization of the form factors, as stated in Appendix A.l. In 
the plots of the branching ratio and FB asymmetries we have used the form factors given 
by Kriiger and Melikov [16]. We have also presented some of the results which one can 
obtain from the form factors given by Dincer and Sehgal [19]. Dincer and Sehgal have used 
a universal form of the form factors as given in Appendix A. 2. The form factors given by 
Kriiger and Melikov [16] and Dincer and Sehgal [19] have one common feature, which differs 
from the earlier formulations of the form factors. The common feature being that away from 
the charmonium resonances the FB asymmetry vanishes. The zero of the FB asymmetries 
in B s — > £ + £~7is a typical feature of the SM, and can be predicted with small theoretical 
uncertainities. This phenomena is similar to the case of the B — > K*£ + £~ decay where the 
zero of the FB asymmetry can be used as a tool for the measurement of Cy /Re(Cg ) [26]. 
Our results also reaffirms the behaviour that FB asymmetry does have a zero in B s — > £ + £~ r y. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



s s 



Figure 2: The differential decay rate (left) and FB asymmetry (right) for B s — > fi + fi 7. 



In our results we have plotted the branching ratio and FB asymmetries as a function of the 
dilepton invariant mass. In the set of graphs, following these, we have also shown the results 
of varying the integrated branching ratio and FB asymmetries against the parameters of the 
2HDM and mSUGRA. For calculating the branching ratio and averaged FB asymmetries in 
the case where \x is the final state lepton we have excluded the region between 0.33 < s < 
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0.55. This region corresponds to the J/ip and ip' resonances. Essentially, in the case of the 
final state lepton being // we have divided the region of the dilepton invariant mass in to two 
regions, defined as; 



Region - I 
Region - II 



4m 2 . < s < 0.33 
0.55 < s < 1-5 



In our analysis we have used a hard photon in the B s 



"7decay and have imposed 



a cut on the photon energy 2 (E 7 ). Where this photon energy cut can be related to the 
parameter 5 by the relation 5 = 2E 1 /rriB- Our estimates have used a photon energy cut 
of Ery = 20MeV. In the case of the final state lepton being r we have used the range 
0.55 <s<l-5. 

The plots of the variation of the differential decay rate and FB asymmetry as a function 
of the dilepton invariant mass (for B s — > /i + /i - 7) are shown in figure 2. Similarly the same 
plot for t in the final state is given in figure 3. As we can see from both figures there can be 
a substantial variation in both these observables when we have an enhanced bsZ coupling. 
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Figure 3: The differential decay rate (left) and FB asymmetry (right) for B s — > t + t 
the 2HDM type-II model. The model parameters are vnn± = 400 GeV and tanfi = 80. 



7 in 



We have presented the results of the branching ratio and averaged FB asymmetry for 
B s — > /i + / u _ 7in Table 1. In the table we have shown our predictions for the SM and with the 
Cio value given by eqn.(23). We have also quoted the results using form factor definitions 
for both Kriiger and Melikov as well as Dincer and Sehgal. As evident from the values given 
in Table 1 the D & S form factors leds to higher values for the branching ratio and averaged 

2 When the photon is "soft" one has to consider both the processes B s — > £ + £~&nd B s — > ^ + ^~7together. 
By taking these two processes together the infrared divergences of B s — > £ + ^~7are cancelled by order a 
corrections of B s — * £ + £~ 



11 



Model 


Observable 


Region - I 


Region - II 


Total 


SM (K & M) 


Branching ratio xlO 9 


0.715 


1.08 


1.796 




< A FB > 


0.265 


0.30 


0.286 


bsZ (K & M) 


Branching ratio xlO 9 


1.93 


4.75 


6.68 




< A FB > 


- 0.05 


-0.07 


-0.068 


SM (D & S) 


Branching ratio xlO 9 


1.2 


1.74 


2.94 




< A FB > 


0.3 


0.45 


0.39 


bsZ (D & S) 


Branching ratio xlO 9 


3.26 


7.21 


10.47 




< A FB > 


-0.071 


-0.11 


-0.098 



Table 1: Integrated decay rate and averaged FB asymmetries for B s — > fi + fi~'j. "K & M" 
refers to the results obtained using the form factors of Kriiger and Melikov [16] whilst "D 
& S" the results obtained using Dincer and SehgaVs form factors [19]. bsZ refers to results 
obtained using the value of Cw as given in eqn.(23). 



FB asymmetries as compared to K & M form factors. We can also observe that with the 
introduction of an enhanced bsZ coupling we have; 



We also observe that; 



Br(B s -> /i + /i- 7 ) 
Br(B s -> h + h~i)sm 



Br(B s -> /i + /i- 7 ) 
Br(B s -> ii + ii-i)sm 

Br(B s — > ^ + /i~7) 
Br(B s -> /i+/i"7)5M 



3.5. 



4m2<s<0.33 



0.55<s<(l-5) 



2.7, 



4.2, 



which means that the enhancement in the branching ratio is more in the high dilepton 
invariant mass (s) region-II. In the case of the FB asymmetries in B s — > yU + /i~7not only 
does the magnitude decrease but there is also a change in the sign of the averaged FB 
asymmetries. 

In table 2 we have quoted the branching ratio and averaged FB asymmetry values for the 
B s — > T + r~7decay in the two models we have considered. We can observe that the enhanced 
Cio Wilson gives an enhancement in the branching ratio of B s — > r + T~7by a factor of 5. 
We can have far greater enhancements in the branching ratio in the 2HDM and mSUGRA 
models. But in all these cases the FB asymmetry tends to decrease. 

In figure 4 we have plotted the integrated branching ratio and averaged FB asymmetries 
of B s — > /i + /!~7in the 2HDM type-II as a function of tanf3 for various values of the charged 
Higgs mass. It can be observed from the figures that the branching ratio tends to increase 
with tan[3 whereas the FB asymmetry tends to decrease in magnitude. We have shown 
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Figure 4: Branching ratio (left) and averaged FB asymmetry (right) for the B s — > ji + ji 7 
decay in the 2HDM. 



Model 


Branching ratio xlO 9 


< A FB > xl(T 2 


SM 


15 


1.6 


bsZ 


72.81 


6 


2HDM 


118.9 


0.33 


mSUGRA 


94 


0.3 



Table 2: Integrated decay rate and averaged FB asymmetries for B s — > t + t~^. In the 
above results we have used the K & M form factors [16]. The parameters for the 2HDM 
model are: vnn± = AOOGeV and tanf3 = 80. Parameters for the mSUGRA model are: 
m = M = 500 GeV, A = 0, tan[3 = 45 and sgn(ii) is taken to be positive. 

similar plots for B s — > r + r~7in figure 5 for the 2HDM type-II. As can be seen from this 
figure there can be an increase in the branching ratio by more than one order of magnitude 
as compared with its SM value. 

In our final set of graphs, figure 6, we have plotted the integrated branching ratio and 
averaged FB asymmetries of B s — > T + r~7as a function of tanfl for various values of the 
unified scalar masses in the mSUGRA model. Here also we have a behaviour similar to 
the 2HDM model where the branching ratio increases with tan(3 but the FB asymmetry 
decreases in magnitude as compared to their respective SM values. 

Radiative decays within the SM have branching ratios ~ 10~ 9 (for /i) and ~ 10~ 8 (for 
r) which can be enhanced by an order in magnitude if some of the new physics comes in 
to play. The enhanced bsZ vertex not only predicts an increase in the magnitude of the 
branching ratio but also predicts a sign change in the averaged FB asymmetry in the case 
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Figure 5: Branching ratio (left) and averaged FB asymmetry (right) for B s — > t + t jin the 
2HDM. 



of B s — ► yU + /i~7relative to the SM. This feature is critically dependent on the nature of 
the form factors used. Therefore future measurements of B s — > £ + £~'-fat hadron collider 
experiments would provide us with more information on the underlying structure of an 
effective Hamiltonian which describes b — > s£ + ^~transitions. 
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A Form Factors 



A.l Kriiger &; Melikov form factors 

The form factors Fy, F TV , Fa and F T a are parameterized as [16]; 

= (24) 

where E 1 is the photon energy. This energy can be related to the dilepton invariant mass 
(in the 5-meson rest frame) as; 

2 \ m z B 
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X 





Figure 6: Branching ratio (left) and FB asymmetry (right) for B s — > r^r^^in the mSUGRA 
model. The mSUGRA parameters we have used are M = 500GeV, A = and sgn(fj) as 
positive. 



0<q 2 <m 2 B , < £ 7 < E n 



For massless leptons the kinematically allowed range would be; 

= m B /2). 

For leptons of mass me, the kinematically allowed range for E 1 would be; 

m\ — Am 2 . 



4m e < q < m B 



, < £ 7 < 



2ms 





F v 


F T y 


F A 


Ft a 


^(GeV- 1 ) 


0.28 


0.30 


0.26 


0.33 


A (GeV) 


0.04 


0.04 


0.30 


0.30 



Table 3: Form factors given by Kriiger & Melikhov [16]. 



A. 2 Dincer & Sehgal form factors 

Dincer and Sehgal have used the universal behaviour for all the form factors given by [19] in 
their definition; 

Fv{Ey) = Fa{E 1 ) = F TV (E y ) = F TA (E,) « — , (25) 

with x 7 = 2E 1 /mB- 
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B Input parameters 



m t = 176GeV , m c = 1.4GeV , m M = 0.105GeV , m T = 1.77GeV , 
m B = 5.26GeV , m b = 4.8GeV , V tb V* = 0.047 , T B = 4.22 x 10~ 13 GeV , 
sin 2 6 w = 0.23 , a = 1/130 , A s = 0.5GeV. , / B = 0.2GeV. 
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